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ABSTRACT 

In the era of high resolution submillimeter interferometers, it will soon be possible to observe the neu- 
tral circumstellar medium directly involved in gas giant planet (GGP) formation at physical scales previously 
unattainable. In order to explore possible signatures of gas giant planet formation via disk instabilities, we have 
combined a 3D, non-local thermodynamic equilibrium (LTE) radiative transfer code with a 3D, finite differ- 
ences hydrodynamical code to model molecular emission lines from the vicinity of a 1.4 Mj self-gravitating 
proto-GGP. Here, we explore the properties of rotational transitions of the commonly observed dense gas tracer, 
HCO + . Our main results are the following: 1. Very high lying HCO + transitions (e.g. HCO + J=7-6) can trace 
dense planet forming clumps around circumstellar disks. Depending on the molecular abundance, the proto- 
GGP may be directly imageable by the Atacama Large Millimeter Array (ALMA). 2. HCO + emission lines 
are heavily self-absorbed through the proto-GGP's dense molecular core. This signature is nearly ubiquitous, 
and only weakly dependent on assumed HCO + abundances. The self-absorption features are most pronounced 
at higher angular resolutions. Dense clumps that are not self-gravitating only show minor self-absorption fea- 
tures. 3. Line temperatures are highest through the proto-GGP at all assumed abundances and inclination 
angles. Conversely, due to self-absorption in the line, the velocity-integrated intensity may not be. High angu- 
lar resolution interferometers such as the Submillimeter Array (SMA) and ALMA may be able to differentiate 
between competing theories of gas giant planet formation. 

Subject headings: planetary systems: protoplanetary disks — planetary systems: formation — line: formation 
— line: profiles — radiative transfer - circumstellar matter 



1. INTRODUCTION 

Since the discovery of 5 1 Pegasus, there have been numer- 
ous detections of Jupiter-sized extrasolar planets through ra- 
dial velocity experiments, transiting of parent stars and di- 
rect imaging (for recent reviews, see Udry, Fischer & Queloz, 
2006, Charbonneau et al. 2006 and Beuzit et al. 2006). 

Concomitant to the problem of characterizing the nature of 
these gas giant planets (GGPs) is developing the theoretical 
constructs which describe the nature of GGP formation. Two 
major theories have been developed concerning the physics of 
GGP formation. Core accretion begins with the formation of 
planetesimals through the collisional coagulation and sticking 
of progressively larger solid bodies in the circumstellar envi- 
ronment. Once ^kilometer sized planetesimals are formed, 
runaway accretion to Mars-sized bodies can occur. When 
the planetary embryos reach a mass of ten Earth masses or 
so, disk gas is accreted dynamically, resulting in GGP for- 
mation. Core accretion is the generally favored mechanism 
for forming Jupiter and Saturn (Pollack et al. 1996; Goldre- 
ich, Li th wick, & Sari, 2004a). The core accretion model has 
been challenged in explaining how a solid core can form on 
timescales less than the disk dissipation times of ~10 -10 7 yr 
(Pollack et al. 1996; Ikoma, Nakazawa, & Emori, 2000). 
However it should be noted that recent models by Rafikov 
(2003) and Goldreich, Lithwick & Sari (2004b) have sug- 
gested a core accretion mechanism that is not incompatible 
with the short-lived disk lifetimes. Additional problems of 
the core-accretion model include resolving theoretical GGP 
core masses with those of Jupiter and Saturn (Mejia, 2004; 
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Saumon & Guillot, 2004). 

As an alternative to core-accretion, the disk instability 
mechanism has been investigated in many models by Boss 
(2001, 2004), Mejia (2004), and Mayor et al. (2005). In 
this theory, marginally unstable gaseous disks contract grav- 
itationally to form GGPs. Models of gravitationally unstable 
disks suggest planets can form on rapid (~10 3 yr) time scales 
(Boss 1997, 1998), but require disks to cool efficiently in or- 
der for gravitationally bound clumps to form. Boss (2004) has 
shown that convective cooling through protoplanetary disks 
can be quite efficient. Hybrid mechanisms for GGP forma- 
tion have been suggested as well (Currie, 2005). 

It is clear that both leading theories in the formation of 
GGPs have their shared successes and challenges in current 
models. Observations of protoplanetary disks will be key in 
constraining the two models. Because the gravitational in- 
stability method involves the accumulation of large clumps 
of cool gas from the circumstellar disk, molecular line obser- 
vations may be helpful in revealing the nature of GGP for- 
mation. Indeed, in the era of high-resolution millimeter and 
submillimeter wave interferometers, clumps of cool and dense 
gas may indeed be directly imageable in nearby circumstellar 
disks. 

Rotational transitions (J+l— > J) in interstellar molecules 
have long been used to better understand the nature of cold gas 
in star-forming environments (for a recent review, see Evans, 
1999). Serving as a proxy for the observationally elusive 
molecular hydrogen (H2), excitation analysis of lines from 
molecules such as CO, CS, HCO + , HCN and others can pro- 
vide diagnostics for the physical conditions in the cold molec- 
ular gas. 

Through the use of submillimeter and millimeter wave in- 
terferometers, direct imaging of cold circumstellar disks have 
recently been made possible (e.g. Patel. et al. 2005, Qi et 
al. 2004). Submillimeter molecular line emission has been 
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FIG . 1 . — Solutions of our radiative transfer codes to test problems "2a" and 
"2b" of van Zadelhoff et al. (2002). The ordinate is fractional populations of 
the J=l and J=4 levels of HCO + and the abscissa is radial distance out from 
center of a 1 -dimensional inside-out collapsing cloud. The solid line is our 
solution, and the dash-dot line is the average solution to the problem taken 
from seven other researcher's codes. 



RH0MAX= -8.3 C0NDIF= 0.3 R= 0.30E+15 




FIG. 2. — Equatorial density contours for hydrodynamic snapshot. Con- 
tours denote changes by a factor of 2 in density, with p max =5.0x 10~ 9 gm 
cnT 3 . The maximum density clump is located at ~4 o'clock in this image. 



used to estimate molecular depletion factors (e.g. Andrews & 
Williams, 2005; van Zadelhoff et al., 2001), as well as pro- 
vide information concerning kinematics in the circumstellar 
environment (e.g. Qi et al. 2003). It is additionally possible 
to constrain disk molecular gas masses through such obser- 
vations (e.g. Hogerheijde et al. 2002). Recent CO (J=6-5) 
observations of TW Hya by Qi et al. (2006) has shown that 
molecular excitation can be quite high due to the dense gas 
present in protoplanetary disks. 

In addition to observational data sets of molecular line 
emission from protoplanetary disks, radiative transfer model- 
ing of the emission patterns can provide powerful constraints 
as to the physical conditions in the molecular gas. As an ex- 
ample, Qi et al. in the aforementioned study of TW Hya used 
a 2D non-LTE Monte Carlo model to derive vertical temper- 
ature distributions in the disk. Similar studies have been per- 
formed to calculate temperature and/or density distributions 
by Kessler (2004), and Semenov et al. (2005), among many 
others. 

In order to fully utilize the new generation of interfer- 
ometers to probe planet formation, submillimeter and mil- 
limeter wave emission modeling of GGP forming disks is 
needed. Progressive works by Wolf & D'Angelo (2005), 
Moro-Martin, Wolf & Malhotra (2005), and Varniere et al. 
(2006) have studied the effects of protoplanets in disks on in- 
frared SEDs through the use of radiative transfer modeling. 
A necessary complementary dataset to these works is models 
of molecular line emission from planet forming circumstellar 
disks. 

In this study, we take the first step in this direction by ap- 
plying a newly developed 3D non-LTE radiative transfer code 
to a model of a gravitationally unstable protoplanetary disk 
(Boss, 2001). We will discuss the emission patterns by way 
of contour maps and emission line profiles. This paper is or- 
ganized as follows: in §2, we discuss the numerical methods 
involved concerning hydrodynamics, radiative transfer, and 
chemistry; in §3 we present synthetic images of HCO + emis- 
sion in the disk; in §4 we discuss emission line profiles; in §5 
we briefly discuss spectral maps and in §6 we summarize. 



2. NUMERICAL METHODS 

2.1. Hydrodynamics 

We have run three dimensional hydrodynamic models using 
a finite-differences code to model the gravitationally unstable 
protoplanetary disk as fully described in Boss (2001). The 
code has been shown to be accurate to second order in space 
and time (Boss & Myhill 1992). The spherical grid is uni- 
formly spaced in the radial direction between 4 and 20 AU. 
The azimuthal grid is uniform, and the polar grid spaced such 
that the resolution grows toward the midplane (with a maxi- 
mum resolution of A#=0.3°). The model includes a central 
protostar of 1 M Q and a disk of 0.091 M Q . The protostar 
wobbles in response to the evolving disk, such that the cen- 
ter of mass of the system is preserved. The initial disk sur- 
face density is assumed to vary with radius as a oc r" 1 / 2 to 
a oc r~ l through the inner disk, and a oc r" 3 / 2 in the outer disk. 
The initial density distribution is seeded with perturbations of 
the form cos(m</>) where m=\, 2,3,4 with an amplitude of 0.01. 
Random noise is included at a lower amplitude. 

2.2. Radiative Transfer 

We build the emergent spectrum by integrating the equa- 
tion of radiative transfer over numerous lines of sight through 
the hydrodynamic snapshot (e.g. Walker, Narayanan & Boss 
1994). The solution to the equation of radiative transfer has 
the numerical form 

I v = J2S V ( Z ) [1 -e~ T » (z) ] e - r " (tot) (1) 

where S v is the source function, and r is the optical depth. If 
local thermodynamic equilibrium (LTE) conditions hold, then 
the source function can be replaced by the Planck function. 
LTE is a fair approximation when the density is much greater 
than the critical density and collisions thermalize the gas, e.g. 
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FIG. 3. — The effect of transition choice on the image: HCO + transitions J=7-6, 4-3, 1-0. The proto-GGP is evident only in the highest HCO + transitions 
due to the high critical density needed to image the dense object. The synthesized images are the mirror image of the density contours in Figure|2| and thus the 
proto-GGP is at ^8 o'clock in the left most panel. The intensity is in units of K-km s -1 and is on a fixed scale for the entire figure. 



where A is the Einstein rate coefficient for spontaneous emis- 
sion, C are the collisional rate coefficients, and u,l correspond 
to the upper and lower states of a given transition. However, 
when considering either low density environments, or molec- 
ular species with high dipole moments (and thus high Einstein 
A coefficients), the approximation of LTE may no longer hold. 
In these situations, the source function must be calculated ex- 
plicitly. 

In order to find the source function, an iterative procedure 
can be employed. We have developed a three-dimensional 
non-LTE radiative transfer code based on the Monte Carlo 
method. The first work to detail a Monte Carlo algorithm for 
non-LTE line transfer was by Bernes (1979). Future work- 
ers improved the algorithms and expanded to two and three 
dimensions (Choi et al. 1995, Juvela et al. 1997 Park & 
Hong 1998, Hogerheijde & van der Tak 2000, Schoier, 2000). 
Monte Carlo techniques in radiative transfer are powerful in 
that they offer a large amount of flexibility for different ge- 
ometries and are easily parallelized. Due to the statistical na- 
ture of Monte Carlo, the solutions include shot noise, and can 
be slow to converge in cases of high (r >100) optical depths. 

The non-LTE radiative transfer code we have developed 
operates on the following principles, adapted from Bernes 
(1979): the goal is to solve for the steady state distribution 
of energy states among the molecules. Once these level pop- 
ulations are known, the source function for a given transition 
can be calculated by: 



tluAu 



(n,B lu -n u B u i) 



(3) 



where B are the Einstein rate coefficients for absorption and 
stimulated emission. However, the problem is circular: the 
level populations in any grid cell depend on the mean intensity 
field through that point. 
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(4) 



The intensity field depends on the emission from other cells 
which in turn is given by their source functions. Hence, it is 
necessary to guess the level populations, solve for the mean 



intensity field, calculate updated level populations and iterate 
until the level populations have converged. 

The radiation field is modeled by photon 'packets' that rep- 
resent many real photons. The number of photons each model 
packet represents is proportional to the Einstein-A for the tran- 
sition and the number of molecules or atoms in the upper state 
of the transition in the emitting grid cell. The photons are 
isotropically emitted in a spontaneous manner with a line fre- 
quency drawn from the line profile function: 



i 



exp 



v— vq — v • n- 



(5) 



where v is the frequency of the emitted photon, vq is the rest 
frequency of the transition, v is the velocity of the emitting 
clump of gas, and a is the standard deviation of the profile 
function. The standard deviation is the Doppler width deter- 
mined by the local kinetic temperature and microturbulent ve- 
locity: 

~2kT 

+ V turb 



a = 



(6) 



We assume a constant microturbulent velocity of 0.2 km s -1 
. The photon then takes a step of a given distance s, passing 
through gas with opacity: 

Ov(dllSt) = KVPdust (7) 



a"' (gas) = -^-4>(v)(niBi u -n u B ul ) 
Air 



(8) 



After passing through this grid cell, the number of real pho- 
tons the model photon packet represents is diminished by a 
factor e~ T due to absorptions. The model photon continues to 
take steps in the same direction until it either leaves the grid 
or the number of photons it represents has become negligible. 
When all of the photons have been emitted, the mean intensity 
is known through each grid point. 

The updated level populations are then calculated with the 
equations of statistical equilibrium 



n, 
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^ n k A kl + ^2 nk(B k iJv + Qz) 

k>l kjl 

which are solved through standard matrix inversion methods. 
With the new level populations in hand, the process can be 
repeated with a new calculated radiation field. These radia- 
tive transfer calculations are then iterated over until the level 
populations are found to converge. Convergence in non-LTE 
radiative transfer simulations depends both on the number of 
model photons realized in the iteration, as well as the optical 
depth. 

The level population calculations are sensitively depen- 
dent on the accuracy of the rate coefficients. We have ob- 
tained our coefficients from the Leiden Atomic and Molecular 
Database 4 (Schoier et al. 2005). In order to independently 
test our radiative transfer codes, we have run the test problem 
of an inside-out collapsing sphere published by van Zadelhoff 
et al. (2002). We present the results of this test in Figure ^ 
and direct the reader to van Zadelhoff et al. (2002) for details 
on the test problem and solutions. 

2.3. The Model and Parameters 

We have run our non-LTE radiative transfer code for HCO + 
rotational transitions through model HR of Boss (2001). In 
this model, a 1 .4 Jupiter-mass dense clump of gas is formed 
through gravitational instabilities in the circumstellar disk. 
The maximum density through the proto-GGP is 5.0 x 10~ 9 g 
cm" 3 , and it orbits at a semimajor axis of ~10 AU. The tem- 
perature through the proto-GGP ranges from ^100-150 K. 

For the radiative transfer calculations, we considered HCO + 
molecular line emission, as well as thermal continuum emis- 
sion with opacities given by Boss & Yorke (1990). The ra- 
diative transfer calculations were run on one snapshot of the 
hydrodynamic model for which we show the midplane den- 
sity contours in Figure |2] The boundary conditions for the 
non-LTE transfer included photons from the 2.73 K cosmic 
microwave background. However, when convolving our re- 
sultant images to beams larger than the grid, we assumed a 
vacuum boundary. Our grid for the calculations was spherical 
in nature with (JV,.,A^,./V e )=(100,512,43). We emitted roughly 
3.11xl0 8 photons per iteration. The calculations typically 
took 2 weeks on 14, 2GHz, AMD-64 processors. 

2.4. Chemistry of HCO + 

Abundances in protoplanetary disks can be wildly variable 
due to the wide range of temperatures and densities involved, 
in addition to variable X-ray and UV fluxes. Our purpose is 
to investigate gross emission features from dense gas clumps 
in circumstellar disks and a full model of the complex chem- 
ical reaction networks in disk environments (e.g. Aikawa & 
Herbst, 1999; Semenov, Wiebe, & Henning 2004) is beyond 
the scope of this paper. We approximate the impact of disk 
chemistry by running models at different HCO + abundances 
and analyzing their imprint on emission features. 

There are commonly two methods of introducing chemi- 
cal depletion in disks: a uniform depletion factor across the 
entire disk, and a jump depletion where grain depletion only 
comes into play below a certain temperature threshold. The 
formation of HCO + depends directly on the abundance of CO, 
and thus the HCO + abundance is assumed to follow the CO 
characteristic depletion. CO freeze-out typically occurs at 
T<22K, thus implying a necessary HCO + depletion at low 

4 http://www.strw.leidenuniv.nl/~moldata 
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FIG. 4. — The effect of angular resolution on the image: HCO + transi- 
tions J=7-6. The left column is at model resolution. The right column is a 
simulated ALMA image at 0.007 " (the most extended baseline of ALMA, 
Bastian, 2002) for a source distance of 140 pc. The synthesized images are 
the mirror image of the density contours in Figure|5] and thus the proto-GGP 
is at ~8 o'clock in each panel. The intensity is in units of K-km s and is on 
a fixed scale for the entire figure. 

temperatures. However, the abundance of HCO + can be al- 
tered from typical interstellar values for a variety of other 
factors: photodissociation can occur in the hot inner regions 
of the disk due to the enhanced UV flux of a young star, or 
through interstellar cosmic rays. We therefore utilize a uni- 
form depletion from standard ISM abundances taken from 
Lee, Bettens & Herbst, (1996). This assumption may limit 
the predicted detectability of massive gas clumps in a circum- 
stellar disk. In a study of the fractional ionization in disks, 
Semenov et al. (2004) find the dominant ion in the interme- 
diate layer and midplane to be HCO + at the radii where the 
densest gas clumps in our models form (10 AU), whereas the 
fractional HCO + drops in the surface layers. By uniformly de- 
pleting the HCO + abundance, we may be decreasing emergent 
flux. 

van Zadelhoff et al. (2001) estimate depletions ranging 
from 0.1-0.001 from standard ISM abundances for disks TW 
Hya and LkCa 15. Dutrey et al. (1997) find an average HCO + 
abundance of 7.4x 10~ 10 /H2, corresponding to a depletion of 
~0.15. We have run three models corresponding to depletion 
factors of 0.5, 0.1 and 0.01. While simple depletion models 
are likely not valid for th e ex tremely dense cores of proto- 
GGPs, as we discuss in § 14.11 the specifics of the chemistry 
in the densest regions of the disk may not contribute signifi- 
cantly to the emergent HCO + flux. 

3. IMAGES 

High angular resolution interferometers such as ALMA 
will be able to achieve unprecedented imaging capabilities of 
Galactic circumstellar disks. Utilizing our radiative transfer 
codes, we have created synthetic intensity contour maps of 
the HCO + emission from the gravitationally unstable proto- 
planetary disk in our models, and present them in this section. 
The parameter space we explore includes rotational transition, 
inclination angle, and abundance. 

3.1. Molecular Transition 

In Figure [5] we show model contours of spectral line inten- 
sity of the protoplanetary disk as viewed through the HCO + 
J=l-0, 4-3 and 7-6 transitions at 1/2 ISM abundance. The 
disk is face on (/=90°) in each image. 

The GGP in our simulation achieves a maximum density 
of n=3.0x 10 15 cm" 3 . At these high densities, and warm tem- 
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FIG. 5. — The effect of inclination on the image: HCO + transition J=7-6. The inclinations are 90° (face on), 45° and 8°. 8° is roughly the minimum inclination 
angle at which the proto-GGP did not get lost in the emission of the disk. The synthesized images are the mirror image of the density contours in FigureEl and 
thus the proto-GGP is at ~8 o'clock in each panel. 
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FIG. 6.— The effect of HCO + abundance on the image. Abundances (from left) are Ixl0" 9 /H 2 , 5xlO _10 /H 2 , 5xlO~ n M 2 . The proto-GGP is only visible in 
the image at highest HCO + abundance. The synthesized images are the mirror image of the density contours in Figure|2] and thus the proto-GGP is at ~8 o'clock 
in the leftmost panel. 



peratures (~100 K), collisions will ensure that most of the 
HCO + molecules are excited well above the ground state. In- 
deed, as our simulations show, the GGP is only fully visible 
at the highest HCO + transitions. While lower transitions in 
HCO + may be able to detect dense clumps in the circumstel- 
lar disk, the densest clumps that may be self-gravitating can 
be identified using high-density tracers. For the HCO + J=7-6 
transition (n cr j t ~10 7 cm" 3 ), the proto-GGP emits quite prodi- 
giously while the clumps in the remaining parts of the disk 
begin to fade. The high J levels of dense gas tracers serve as 
an efficient method of filtering out low density gas that may 
not be directly associated with the protoplanet. We have not 
explored any transitions beyond HCO + J=7-6. It may be, how- 
ever, that the proto-GGP is even more distinct at higher (THz) 
HCO + transitions. For the transitions we have modeled, the 
densest clumps in the circumstellar disk are most visible in 



HCO + J=7-6. We will explore the effects of inclination and 
abundance primarily in this transition. 

With its most extended baseline, ALMA will be able to 
achieve a spatial resolution of 0.007" at HCO + (J=7-6, v =624 
GHz) (Bastian 2002). At the distance of Taurus Molecular 
Cloud, ~140pc, this angular resolution (1 AU) is quite com- 
parable to our model resolution of ^0.5 AU. In Figure 0] we 
have simulated an observation of our circumstellar disk by 
setting it at a distance of 140 pc and convolving it with the 
0.007" ALMA beam. We assume no lost flux and a circular 
Gaussian beam. The proto-GGP and other dense clumps are 
quite visible. As an example, in one hour of integration, the 
64-element ALMA array should be able to image the dense 
protoplanet in Figure |4] with a signal to noise ratio of ~3- 
4. Additionally, as we will show in § signatures of dense 
clumps may be evident even in single-dish sub-mm telescopes 
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FIG. 7. — HCO + J=7-6 spectral line profile at model resolution through 
the proto-GGP with optical depth overplotted. The central density is 
~10 I5 cirT 3 , and the line is completely self-absorbed at line center. The 
emission is primarily due to radiatively pumped gas in the outer layer of the 
proto-GGP. 




FIG. 9. — HCO + J=7-6 spectral line profile of proto-GGP, convolved to the 
beam of different telescopes. The assumed abundance is 1 X 10~'/H2 (0.5 de- 
pletion from standard ISM values, Lee, Bettens & Herbst, 1996). In all cases, 
the simulated telescope is pointed directly at the proto-GGP. The telescopes 
simulated are ALMA, the SMA and a 10m class single dish. The correspond- 
ing optical depth, t, is overplotted, with values on the right side of each box. 
The disk is face-on. 
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FIG. 8. — HCO + J=7-6 spectral line profile at model resolution through the 
dense clump at position x = 5 AU, y = -3 AU (for reference with Figure 151. 
The central density of the clump is ~10 15 citT 3 . Some emission from the 
central regions escapes, and the line is not completely self-absorbed. 

through signatures in the emission line profile. 

3.2. Inclination 

As we will discuss in § 14.11 the emission from the proto- 
GGP is largely from the outer layers. Consequently, the in- 
clination angle of the disk does not affect the results signifi- 
cantly. As seen in Figure [5] even at nearly edge-on inclina- 
tions, the proto-GGP is still visible at HCO + J=7-6. However, 
due to increasing column through the disk in low-inclination 
scenarios, other dense clumps in the disk begin to emit at sim- 
ilar intensities as the proto-GGP, thus confusing observations. 
Because of this effect, we have found that the minimum angle 
in our models that the proto-GGP is visible is i ~8°. Inclina- 
tion effects can quickly wash out signatures of the proto-GGP 
in the spectral line profile (see |4}. 

3.3. Abundances 



We have modeled the HC O + abundances as uniform, mod- 
ulo depletion factors (§ 12.41 1. In Figure|5] we present face-on 
images of the protoplanetary disk at 0.5, 0.1 and 0.01 deple- 
tion factors from standard ISM abundances. The proto-GGP 
is no longer visible at lower abundances. This holds for all in- 
clinations, as well as all lower HCO + transitions. With less 
emitting molecules along the line of sight, the proto-GGP 
emission is no longer able to dominate over the other dense 
clumps in the disk. However, the proto-GGP does, in fact, 
emit similar peak line temperatures as other dense clumps 
in the cloud. Why, then, does the protoplanet not appear to 
emit brightly in the image? The image shows velocity inte- 
grated line intensity - the total area under the emission spec- 
trum. Much of the emission from the proto-GGP is self ab- 
sorbed, and the majority of the emission comes from radia- 
tively pump ed gas in the outer layers at low velocity disper- 
sion (§ 14. It . This effect reduces the velocity integrated inten- 
sity. 

4. LINE PROFILES 

4.1. Non-LTE Effects: Radiative Pumping in the Vicinity of 
Dense Gas Clumps 

In Figure we have plotted the HCO + (J=7-6) emission 
spectrum through the proto-GGP in the model with 0.5 de- 
pletion factor. The large optical depths cause the line flux to 
be ^0 K at line center. The emission from the proto-GGP 
does not originate at the core of the dense clump of gas, but 
rather in the more diffuse outer layers. While the gas in the 
dense core of the proto-GGP is in LTE, the emission from this 
gas suffers heavy extinction owing to the high optical depths 
at line center. However, the density through the proto-GGP 
drops off quickly with radius allowing gas in the outer layers 
(where r <1) to be radiatively pumped by ^1500 K gas near 
the protostar. The strong emission line temperatures, then, 
originate in radiative excitation of lower density non-LTE gas, 
rather than from the core of the cold proto-GGP, and conse- 
quently does not reflect the kinetic temperature of the emitting 
gas. 

The emission pattern of heavy self-absorption at line cen- 
ter, and radiatively pumped gas at the surface is characteristic 
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FIG. 10. — HCO + J=7-6 spectral line profile of proto-GGP, convolved to 
the beam of different telescopes. The assumed abundance is 5xl(r 10 /H2 
(0.1 depletion from standard ISM values, Lee, Bettens & Herbst, 1996). In 
all cases, the simulated telescope is pointed directly at the proto-GGP. The 
telescopes simulated are ALMA, the SMA and a 10m class single dish. The 
corresponding optical depth, r, is overplotted, with values on the right side 
of each box. The disk is face-on. 

only of the densest clump of gas in our models, e.g. the self- 
gravitating proto-GGP with central density ~10 15 cm" 3 . As 
it is these densest clumps that serve as antecedents to GGPs, 
the line profile of heavy self absorption at line center com- 
bined with offset emission may be characteristic of GGPs in 
formation. The optical depth in other dense clumps in the cir- 
cumstellar disk is sufficient to produce self-absorption in the 
line profiles, but not to the same degree as seen toward the 
proto-GGP. As an example, in Figure [8] we plot the emer- 
gent HCO + J=7-6 spectra at model resolution and depletion 
factor 0.5 from a clump in the face-on disk (located at x = 
5 AU, y = -3 AU for reference with Figure |6) with central 
density ~ 1.75x10 cm" 3 . The optical depth at line center 
only reaches r»l near the dense core of the clump, allow- 
ing significant emission from both the LTE core, as well as 
the radiatively pumped outer layers of the dense clump to es- 
cape. In contrast, the optical depth at line center at the core 
of the proto-GGP reaches a total value of several thousand, 
rendering the majority of the proto-GGP optically thick. 

The nature of emission from the proto-GGP has impli- 
cations concerning the assumed chemistry in these models. 
High optical depths through the densest regions of the proto- 
GGP prevents radiation from emerging. For example, even in 
the model with lowest optical depths (0.01 depletion factor), 
only ~5% of the emission at the line peak originating from 
the center of the proto-GGP reaches the observer. This pales 
in comparison to the emission from the radiatively pumped 
gas in the outer layers of the proto-GGP which has an order 
of magnitude higher source function than the LTE gas at the 
core. The emission due to radiative pumping of the less dense 
gas in the outer layers of the proto-GGP will tend to mask 
more complex chemistry that may be occurring in the denser 
regions. 

4.2. Effects of Abundance and Resolution 

In order to further quantify the effect dense clumps have on 
the emission line profiles, we have convolved the model re- 
sults with circular Gaussian telescope beams of different sizes 
(corresponding to the diameter of the telescope and the fre- 
quency of interest) and plotted the synthesized emission line 
spectra for three different HCO + abundances. 




FIG. 1 1 . — HCO + J=7-6 spectral line profile of proto-GGP, convolved to 
the beam of different telescopes. The assumed abundance is 5xlO~"/H2 
(0.01 depletion from standard ISM values, Lee, Bettens & Herbst, 1996). In 
all cases, the simulated telescope is pointed directly at the proto-GGP. The 
telescopes simulated are ALMA, the SMA and the 10m class single dish. 
The corresponding optical depth, r, is overplotted, with values on the right 
side of each box. The disk is face-on. 



In Figure[5]E3and^2 we ^ ave pl° ttec l tne model emission 
line profile as viewed by ALMA, the SMA and 10m class 
single-dish telescopes [e.g. Caltech Submillimeter Observa- 
tory (CSO) or the Heinrich Hertz Submillimeter Telescope 
(SMT)]. With each spectral line, we have overplotted the op- 
tical depth (r) as a function of velocity. The model disks 
are viewed face on. The line profiles are found to change 
significantly as a function of abundance. However, in each 
instance the line profiles through the proto-GGP are heavily 
self-absorbed. 

In the model with a 0.5 depletion factor, the emission line 
through the proto-GGP, both at model resolution and that of 
ALMA, is self-absorbed to zero flux at line center. This is ev- 
ident from the location of the peak in the optical depth profile. 
The asymmetry in the line profile (~l-2 km s" 1 of line center) 
is due to motion along the line of sight of the emitting atmo- 
sphere. As the beam size gets larger, more emission from the 
disk is included, and the effects of self-absorption begin to be 
washed out. The model with 0. 1 depletion factor is qualita- 
tively very similar to the 0.5 depletion factor case. There are 
minor differences involving the emission of different veloc- 
ity components in the proto-GGP, but generally the results are 
the same. It is notable, however, that the peak line tempera- 
ture has dropped by a factor of 2-3. Self-absorption features 
become more pronounced in the model with 0.01 depletion 
where an inverse p-cygni profile is observed. In this lowest 
abundance model, the column of emitting molecules is low 
enough that the emission can no longer fill in the absorption 
trough. The absorption is washed out for angular resolutions 
>0.1". The emission from the rest of the disk, once folded 
into the beam, quickly counteracts the self-absorption feature 
and a single peaked Gaussian line is observed. 

At all abundances in our models, the self absorption through 
the proto-GGP reduces the flux at line center to or negative 
values. Conversely, non self-gravitating dense clumps in our 
models only show moderate self-absorption. It may be that 
complete self-absorption to zero or negative flux is a signa- 
ture of GGP formation discernible by ALMA. As the disk 
becomes more inclined, the self-absorption features become 
less discernible as the emission line broadens. 
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FIG. 12. — Spectral map of model with 1 X 10 '/Lb abundance. The disk is face on, and the beam is 0.007", corresponding to the longest ALMA baseline. The 
x-axis in each plot is offset velocity (km s ) and the y-axis in Rayleigh-Jeans Temperature (K). The title in each panel refers to the offset position in cm. The 
map is made from the mirror image of the density contours presented in Figurel2l 



5. SPECTRAL MAPS 

In Figure we have taken spectra along various lines of 
sight through the face on protoplanetary disk with 0.5 deple- 
tion, convolved to a 0.007" ALMA beam. The proto-GGP 
stands out in emission quite high above the rest of the disk. 
This was evident as well in the contour emission map (Figure 
®. 

Spectral mosaics are particularly useful for identifying 
GGPs at low HCO + abundance levels (e.g. FigurelT3l where 
the proto-GGP may make itself evident through its bright 
emission line and self-absorbed profile. This type of emis- 
sion is missed in a velocity integrated intensity contour map, 
however, owing to self-absorption in the line. 

6. CONCLUSIONS 

In this paper, we have presented non-LTE radiative transfer 
calculations of HCO + rotational line emission in a gravita- 
tionally unstable protoplanetary disk. Our models suggest the 
following: 

1 . Dense gas clumps associated with proto-GGPs formed 
via disk instabilities may be observable with interfer- 
ometers (e.g. ALMA) in molecular transitions with 
high critical density (e.g. HCO + J=7-6). The emission 
is dominated by radiatively pumped gas in the outer lay- 
ers of dense clumps, and thus does not necessarily re- 
flect the physical temperature of the emitting gas. 



2. The emission lines arising from the densest self- 
gravitating gas clumps in the protoplanetary disk are 
completely self-absorbed to zero or negative flux. The 
line profiles from other dense gas clumps are self ab- 
sorbed, as well, although to a lesser degree. 

3. Proto-GGPs often appear as bright peaks in spectral line 
maps, making them easier to identify in the presence of 
extended disk emission. 

4. ALMA will have the angular resolution and sensitivity 
necessary to directly image proto-GGPs in formation in 
nearby circumstellar disks. 

More complete chemical models and reaction networks will 
be incorporated in our models in future works. 
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